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29  av.  de  la  division  Leclerc,  92322  Chatillon  cedex,  France 
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email  : lionel.larcheveque@onera.fr 

Abstract.  MILES  and  traditional  LES  computations  of  a high  subsonic 
flow  over  a deep  and  open  cavity  using  wall  function  and  2D-3D  domain  cou- 
pling strategy  have  been  carried  out  ( Rei  = 6.7  x 101 * * * 5).  Both  Reynolds  and 
phase  averages  have  been  computed.  Results  show  especially  good  agree- 
ment with  experimental  data. 


1.  Introduction 

High-speed  flow  over  a cavity  generates  a complex  flowfield  whose  dynamics 

are  governed  by  various  physical  mechanisms  such  as  shear  layer  instabili- 

ties, acoustic  forcing,  etc.  This  class  of  flows  is  not  only  of  fundamental  in- 
terest, but  also  occurs  in  many  aerodynamic  configurations,  such  as  weapon 
bay  or  wheel  wells,  where  large  pressure  fluctuations  induced  by  this  flow 
are  known  to  result  in  possible  structural  damages.  The  capability  of  nu- 
merically predicting  these  flows  in  a time-accurate  fashion  is  consequently 
an  important  issue.  In  the  past  twenty  years,  numerous  RANS  computa- 
tions, mostly  supersonic,  have  been  carried  out.  The  increase  of  available 
computational  resources  yet  allows  LES  computations  of  such  configura- 
tions. 

The  studied  cavity  has  a length  L to  depth  D ratio  of  0.42,  and  therefore 
belongs  to  the  deep  cavity  category.  There  is  no  reattachment  of  the  flow 
on  the  floor  (open  cavity).  The  Reynolds  number  (based  on  the  inflow 
velocity  Uoo  and  the  length  L of  the  cavity)  and  the  Mach  number  are 
6.7  x 105  and  0.8  respectively.  This  configuration  has  been  chosen  because 
of  the  availability  of  an  experimental  data  base  compiled  by  Forestier  et 
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al.[l],  which  was  specially  designed  for  the  accurate  evaluation  of  DNS  and 
LES  computations.  These  data  include  averaged  and  phase  averaged  (20 
phases)  velocity  measurements  in  various  locations  along  the  channel  and 
inside  the  cavity.  According  to  this  experience,  the  flow  is  characterized  by 
a fundamental  frequency  of  1975  Hz.  Spectrum  exhibits  multiple  harmonics 
of  decaying  SPL  levels  (see  left  part  of  the  figure  4).  The  length  of  a period  is 
not  constant  (jitter),  the  standard  deviation  being  approximatively  equal  to 
1%  of  the  mean  value.  During  one  period,  shadowgraphs  and  phase  averaged 
velocity  reveal  the  existence  of  three  coherent  vortices  (see  left  part  of  the 
figure  8).  Shadowgraphs  also  show  the  presence  of  strong  pressure  waves. 

2.  Numerical  methods 

Spatial  discretization  is  accomplished  using  a finite  volume  scheme.  Two 
categories  of  computations  have  been  carried  out.  The  preliminary  results 
were  obtained  without  the  use  of  any  subgrid  model,  an  approach  known 
as  MILES,  introduced  by  Boris  et  al.[2],  and  which  is  based  on  the  idea 
that  an  upwind  scheme  is  able  to  mimic  the  dissipative  behaviour  of  the 
small  structures.  According  to  this  technique,  the  convective  fluxes  are  eval- 
uated by  means  of  an  upwinding  formulation  derived  from  the  AUSM+(P) 
scheme  originally  proposed  by  Edwards  and  Liou[3j.  With  the  intention  of 
getting  to  the  root  of  this  study,  a more  traditional  LES  approach  was  also 
tested,  using  the  mixed-scale  model  of  Sagaut  coupled  with  a selective  func- 
tion (Sagaut  and  Troff[4]).  In  that  case,  the  second  order  spatial  scheme  is 
centered  in  order  to  ensure  minimal  numerical  dissipation.  To  prevent  the 
apparition  of  oscillations,  a detector  proposed  by  Mary[5]  locally  switches 
if  necessary  to  the  former  scheme.  Explicit  temporal  integration  is  achieved 
using  a third  order  Runge-Kutta  scheme.  In  order  to  greatly  reduce  the 
computational  cost,  an  instantaneous  log  law  is  used  at  the  wall.  Periodic 
boundary  conditions  are  set  in  the  spanwise  direction.  Velocity  profile  is  im- 
posed at  the  entrance  of  the  channel,  with  additional  white  noise  velocity 
fluctuations.  Both  mean  and  fluctuating  velocities  are  adjusted  in  concor- 
dance with  experimental  measurements.  Non-reflective  outflow  boundary 
conditions  are  used  on  the  left  extremity  of  the  channel. 

The  computational  domain  extends  from  -1L  to  5L  in  the  streamwise 
direction  and  from  -2.4L  to  2L  in  the  vertical  direction,  the  origin  being 
set  in  the  leading  edge  of  the  cavity  (see  figure  2 for  a glimpse  of  the  com- 
plete domain).  The  computation  is  done  on  a multibloc  structured  mesh, 
with  2D-3D  coupling  strategy.  The  upper  side  of  the  channel,  in  which 
the  flow  is  almost  two-dimensional,  is  discretized  in  the  spanwise  direction 
with  less  cells  than  the  lower  part  of  the  channel.  The  remaining  part  of 
the  computational  domain  (including  the  cavity)  is  fully  3D.  Two  meshes 
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have  been  employed,  a coarse  one  with  800,000  cells  and  a fine  one  with 
approximatively  twice  this  number,  obtained  by  multiplying  the  number  of 
cells  in  both  the  streamwise  and  normal  directions  by  1.5,  the  number  of 
cells  in  the  spanwise  direction  remaining  constant.  The  MILES  approach 
was  tested  only  on  the  coarse  mesh.  The  following  table  summarizes  the 
typical  size  of  the  cells  in  wall  unit  for  the  LES  computations. 


Mesh 

streamwise 

spanwise 

-L  to  wall 

Coarse 

150  - 400 

50  - 200 

20-30 

Fine 

100  - 300 

50  - 200 

15  - 20 

Due  to  the  restrictive  CFL  conditions,  dimensionless  time-steps  of  9.1  x 10-4 
on  the  coarse  mesh  and  7.28  x 10~3 4  on  the  fine  one  have  been  adopted. 

As  seen  in  figure  1,  the  procedure  used  to  compute  phase  averages  con- 
sists in  splitting  the  interval  between  five  consecutive  peaks  of  filtered  pres- 
sure into  80  sub-intervals.  The  goal  of  such  an  approach  is  to  provide  against 
dephasing  due  to  the  jitter.  For  all  the  computations  the  standard  devia- 
tions of  the  length  of  the  period  were  found  of  the  same  magnitude  as 
the  experimental  one.  Tests  with  experimental  data  have  proved  than  with 
up  to  four  by  four  periods  splitting,  the  effect  of  the  jitter  is  negligible. 
Therefore  this  upper  value  has  been  adopted.  The  averaging  time  for  both 
Reynolds  and  phase  averages  corresponds  to  13  x 4 = 52  cavity  cycles. 

3.  Numerical  results  and  comparison  with  experimental  data 

Figures  2 and  3,  based  on  instantaneous  data  issued  from  the  fine  mesh, 
show  coherent  structures  highlighted  using  Q criterion  (¥_u  tensor  2nd  in- 
variant) combined  with  a pseudo- Schlieren  view.  The  main  characteristics 
of  the  flow  are  recovered,  with  the  existence  of  three  major  structures  (the 
two  vortical  tubes  on  the  left  of  the  figure  3 are  about  to  pair)  and  pressure 
waves  with  similar  shapes  compared  to  the  experimental  shadowgraphs. 

The  experimental  and  computed  spectra  are  shown  in  figure  4.  The 
concordance  is  really  good,  with  an  accurate  prediction  of  both  the  mean 
value  and  peaks  levels  up  to  26000  Hz  ( 12th  harmonic).  The  frequencies  of 
the  peaks  are  slightly  over-evaluated  (2.5%),  but  it  has  to  be  noticed  that 
the  spectral  resolution  is  about  2.2%  of  the  fundamental  frequency.  This 
spectrum  is  obtained  from  fine  mesh  data.  On  the  coarse  mesh,  the  levels 
are  decaying  faster,  especially  with  the  MILES  computation,  and  only  the 
first  four  peaks  are  present,  with  the  same  levels  as  for  the  fine  mesh. 

The  evolution  of  the  incompressible  momentum  thickness  6 along  the 
shear  layer  is  plotted  in  figure  5.  All  the  three  computed  graphs  are  close 
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of  the  experimental  data,  with  the  exception  of  an  earlier  effect  of  the 
trailing  edge.  MILES  and  LES  computations  on  the  coarse  mesh  give  similar 
results,  with  a slightly  underestimated  growth  rate  in  the  initial  region 
(5  < x/9r  < 20).  Using  the  fine  mesh,  we  obtain  initially  a higher  growth 
rate,  but  the  layer  grows  weakly  in  the  region  40  < x/9r  < 60,  where  the 
flow  is  supposed  to  be  near  to  the  equilibrium.  In  the  first  region,  Forestier 
et  al.[6]  suppose  that  the  fast  growth  may  be  explain  by  an  early  pairing 
of  small  vortical  tubes,  a mechanism  similar  to  the  “collective  interaction” 
introduced  by  Ho  and  Huang[7j.  Therefore  a plausible  explanation  for  the 
higher  growth  rate  may  be  that  because  of  the  refinement  of  the  mesh, 
small  structures  in  the  vicinity  of  the  leading  edge  and  their  dynamics  are 
more  accurately  described. 

Figures  6 shows  the  evolution  of  the  vorticity  thickness  4,  defined  by: 
<f(z)  = {U+  - U~)  / [max*  (§f)] 

where  U+  is  the  velocity  in  the  center  of  the  channel  and  U~  is  the  velocity 
inside  the  cavity,  here  considered  equal  to  0.  Beyond  x/9r  — 25,  computa- 
tions and  experience  are  well  correlated,  especially  when  using  a subgrid 
model.  Here  again,  there  is  an  earlier  influence  of  the  trailing  wall.  Before 
x/9r  — 25,  the  evaluation  of  the  experimental  values  of  4 is  of  a poor  accu- 
racy because  of  the  sharp  vertical  velocity  gradient  resulting  from  the  small 
thickness  of  the  shear  layer  and  the  coarse  size  of  the  measure  mesh. 

Figure  7 presents  Reynolds  and  phase  averaged  streamwise  velocity 
and  2D  turbulent  kinetic  energy  profiles  at  a location  corresponding  to 
x/6r  = 50,  in  the  middle  of  the  self-similarity  growth  region.  Focusing  on 
the  streamwise  velocity,  we  see  that  the  LES  profiles  associated  with  the 
coarse  and  the  fine  meshes  are  almost  indistinguishable.  MILES  profiles 
differ  mainly  inside  the  cavity,  with  smaller  speed.  However  for  each  case 
the  agreement  with  experimental  data  is  good.  If  we  consider  the  turbulent 
kinetic  energy  profiles,  the  two  types  of  LES  calculations  are  still  closely 
connected.  On  Reynolds  average,  MILES  computation  generates  more  tur- 
bulent kinetic  energy.  On  the  contrary,  LES  leads  to  slightly  underestimated 
turbulent  kinetic  energy  peaks,  but  the  shapes  are  correct.  The  two  phases 
selected  show  deficit  in  energy  for  each  case,  when  the  locations  of  the  peaks 
are  well  predicted. 

Figure  8 shows  a comparison  of  the  loci  of  the  coherent  structures  be- 
tween experimental  data  and  numerical  results  on  the  fine  mesh,  for  four 
different  phases.  The  structures  are  highlighted  using  the  Weiss  criterion 
(equivalent  to  2D  Q criterion).  If  we  consider  that  the  experimental  mea- 
surement mesh  is  really  coarse  compared  with  the  computing  mesh  and  that 
this  difference  between  the  two  meshes  has  an  important  effect  on  the  com- 
putation of  Vu,  it  stands  to  reason  that  the  overall  agreement  is  excellent. 
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Each  experimental  structure  can  be  linked  to  a numerical  structure  with 
quite  correct  location  and  size,  except  perhaps  for  the  less  energetic  ones 
(see  structure  3 in  phase  11).  The  results  of  the  computation  on  the  coarse 
mesh  (MILES  and  LES),  which  are  not  shown  here,  are  closely  similar  to 
the  results  on  the  fine  mesh.  As  a remark,  we  can  notice  that  the  figures  2 
and  3 correspond  to  a time  which  takes  place  into  the  first  phase. 

4.  Conclusion 

All  the  computations  presented  here  show  good  agreement  with  the  experi- 
mental data.  This  similarity  despite  various  numerical  methods  and  meshes 
may  be  explain  by  the  energetic  predominance  of  the  low  frequencies  (and 
therefore  the  large  structures).  These  simulations  have  then  proved  the 
(MI)LES  techniques  capability  of  studying  this  category  of  cavity  problems, 
and  therefore  the  possibility  to  use  them  to  improve  our  understanding  of 
the  flow  dynamics  inside  and  above  the  cavity. 
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Figure.  I.  Procedure  to  compute  phase  averages  : localization  of  maxima  Pmax  for 
filtered  pressure  fluctuations  and  splitting  of  the  interval  between  Pmax  and  Pmax  into 
80  sub-intervals  with  equal  duration. 
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Figure  2.  View  of  the  computational  domain.  Structures  are  visualized  using  iso-surface 
of  Q equal  to  10  (^t)  combined  with  a Schlieren-like  picture  in  the  background. 


Figure  3- 


Closeup  view  of  the  mixing  layer  above  the  cavity:  Q value  equal  to  20 
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Figure  4 • Experimental  (left)  and  computed  (right)  spectra.  The  computation  is  done 
on  the  fine  mesh  with  subgrid  model  (LES). 


Figure  5.  Momentum  thickness  6 of  the  mixing  layer  normalized  with  initial  experimental 

momentum  thickness  9r. : MILES  with  coarse  mesh; : LES  with  coarse  mesh; 

: LES  with  fine  mesh;  o;  experimental  data. 


Figure  6.  Vorticity  thickness  <5  of  the  mixing  layer  normalized  with  L. : MILES  with 

coarse  mesh; : LES  with  coarse  mesh;  : LES  with  fine  mesh;  o:  experimental 

data. 
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Figure  7.  Mean  streamwise  velocity  (left)  and  2D  turbulent  kinetic  energy  (right) 
profiles  at  i = 0.6L:  Reynolds  average  (top),  phase  average  for  phase  1 (middle)  and 

phase  average  for  phase  11  (bottom). : MILES  on  coarse  mesh; : LES  on 

coarse  mesh;  : LES  on  fine  mesh;  o:  experimental  data. 
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Figure  8.  Phase  averaged  coherent  structures  isolated  using  Weiss  criterion  (levels  from 
1 to  25  by  increment  steps  of  2,  normalized  using  (^f2-)2):  experimental  data  (left)  and 
LES  on  fine  mesh  (right)  for  phases  1,  6,  11  and  16. 


